Nonequilibrium thermal transport and its relation to linear response 
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We study the real-time dynamics of spin chains driven out of thermal equilibrium by an initial 
temperature gradient Tl ^ Tr. We demonstrate that the nonequilibrium energy current saturates 
fast to a finite value if the linear-response thermal conductivity is infinite, i.e. if the Drude weight D is 
nonzero. Our data suggests that a nonintegrable dimerized chain might support such dissipationless 
transport [D > 0). We show that the steady-state value of the current for arbitrary Tl ^ Tr is 
completely determined by the linear conductance. Inhomogeneous systems exhibiting different bulk 
parameters as well as Luttinger liquid boundary physics induced by single impurities are discussed 
shortly. 
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Introduction — One-dimensional (Id) electronic sys- 
tems are realized in carbon nanotubes and individual 
polymer molecules and provide an approximate descrip- 
tion of strongly anisotropic 3d materials. It has been 
known for many years that Id systems can support un- 
usual correlated electron phenomena such as Luttinger 
liquid physics. However, electrical and thermal transport 
in real materials are usually not governed by the free low- 
energy Luttinger liquid fixed point but by an interplay 
between dangerously irrelevant operators scattering the 
currents and conserved quantities protecting them [ij-yj. 

In order to connect to actual experiments, it is thus es- 
sential to study generic microscopic models. Over the last 
decades a significant number of works [1 - 11 1 investigated 
equilibrium charge (or spin) transport properties. In par- 
ticular, the question whether or not so-called integrablc 
models, which possess a complete set of local conserved 
quantities, can support dissipationless currents at finite 
temperature was addressed extensively. Less is known 
about the quantitative effects of integrability-breaking 
perturbations which are naturally present in any experi- 
mental system, and even the qualitative question whether 
the linear-response conductivity of a nonintegrable model 
can still be infinite is not answered conclusively [l2j ]. 
While experimental measurements of thermal transport 
driven by a temperature gradient in quasi-ld spin sys- 



tems already exist 1 
this theoretically [ij 



14L only a few works investigate 
15|420(. Studying nonequilibrium 
thermal (or charge) transport is complicated in general - 
one reason being that is not even clear whether the long- 
time dynamics can be described by a low-energy theory - 
and constitutes one of the most active areas of research in 
strongly correlated condensed matter physics [21TI29I l32| . 

The primary goal of our work is to obtain quantita- 
tive results on steady-state energy flow both near and 
far from equilibrium and to understand the effects of in- 
tegrability and correlations. This is motivated by the ex- 
periments listed above and by recent technical advances 
in dynamical simulations [111 ]. As a prototypical model 
we consider a XXZ spin- 1/2 chain in presence of two per- 



turbations (dimerization and a staggered magnetic field) 



311. We demonstrate that 



which break integrability [3C 
the nonequilibrium energy current driven by a tempera- 
ture gradient Tt, ^ Tr relaxes fast to a finite steady-state 
value if the linear-response thermal conductivity is infi- 
nite, i.e. if the Drude weight D is nonzero. Our data 
indicates that the dimerized chain might support such 
dissipationless transport (D > 0) despite the fact that it 
is nonintegrable (D can be extracted from the asymptote 
of the equilibrium energy current correlation function , 
and we cannot exclude that the latter decays on a hidden 
large time scale) . We show that for a large class of prob- 
lems the steady-state current takes the functional form 
/(Tl) — /(Tr). This implies that nonequilibrium ther- 
mal transport is entirely determined by linear response 
- / can simply be obtained by integration of the equi- 
librium conductance &rf- We demonstrate that at low 
temperatures the gapless integrable XXZ chain as well 
as the quantum Ising model exhibit universal nonequi- 
librium behavior conjectured by conformal field theory 
32 . 33|. Finally, we study inhomogeneous systems fea- 
turing different bulk interactions as well as Luttinger liq- 
uid physics [3~ij induced by an impurity at the interface. 

Thermal non- equilibrium setup - We aim at inves- 
tigating the real-time dynamics of the energy current 
(JE(n,t)) through a one-dimensional infinite lattice sys- 
tem driven out of equilibrium by an initial sharp tem- 
perature gradient Tr, ^ Tr. Our main focus is to study 
the long-time behavior of (Js(n, t)} and specifically the 
question how it relates to linear-response thermal trans- 
port properties. As a prototypical model, we consider 
a chain of interacting spin-1/2 degrees of freedom S^' y ' z 
governed by local Hamiltonians 
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or equivalently spinless Fermions through a Jordan- 
Wigner transformation. By choosing the couplings J n , 
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FIG. 1. Energy current flowing between two semi-infinite 
spin- 1/2 chains which are initially in thermal equilibrium at 
different temperatures Tl,r and coupled at time t = and 
position n = 0. (a) Integrable XXZ chain with z-anisotropy 
A. The behavior in the gapped phase A > 1 is similar, (b) 
Nonintegrable dimerized XXZ chain where the coupling on 
every second bond is reduced by A. The system is gapped, 
(c) XXZ chain in presence ol a staggered field b rendering 
the model nonintegrable. A gap opens for b » 0.3. Despite 
the fact that the local energy density h(n,t) does not relax 
(Inset), the current saturates fast to a unique finite value ex- 
cept for b > 0. We attribute this to a finite linear-response 
thermal Drude weight of both the pure XXZ chain and the 
nonintegrable dimerized chain (see Figure • This indicates 
an intricate relation between nonequilibrium and linear ther- 
mal transport properties. 



A„, and b n appropriately: 
1 n odd 



J' n 



A n even 



A„ = A , b n = 



(-!)"& 
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we can study systems which are gapless or gapped and 
- as a key aspect of this work - investigate the role of 
integrability. For A = 1 and 6 = 0, Equation ([T]) can be 
diagonalized via Bethe ansatz [3o| ; the model is nonin- 
tegrable otherwise. The spectrum is gapless for |A| < 1 
and gapped for A > 1. A gap opens for A < 1 or suffi- 
ciently large b > b c ; b c > if -1 < A < -l/y/2 [3l|,|36|. 

Thermal nonequilibrium is introduced via the following 
protocol: We initially consider two seperate semi-infinite 
chains (N -> oo), H = H L + H R = J2n=-N/2+i h n + 
'Y^a 1 h n , each being in thermal (grand-canonical) 
equilibrium at temperatures Tl and Tr. The corre- 
sponding density matrix factorizes, po = pl ® Pr., Pi = 
exp(- J H i /T i )/Trexp(-ir i /T i ). At time t = 0, the chains 
are coupled through ho, and the time evolution of p Q is 
computed w.r.t. H = Hq + ho- The energy current is 



FIG. 2. Linear-response energy current correlation func- 
tion whose long-time asymptote determines the Drude weight 
through Equation (J5]). For the integrable XXZ chain (A = 1, 
6 = 0), the global energy current Je — Yl n ^ E ( n ) ^ s con- 
served; thus, (Je^Je) = (Je(0)Je), and DMRG can be 
compared with exact Bethe ansatz results (symbols). (a,b) 
Nonintegrable dimerized chain. The current correlation func- 
tion seems to saturate at a finite value (or decays on a hidden 
large time scale), indicating a finite Drude weight D > 0. (c) 
The data in presence of a staggered field b > is consistent 
with D = 0. 



defined by a continuity equation [8[ , 

d t h n = Je(ji) — J E (n+l) => Jb(ji) = h„] , (3) 

and its time evolution is simply given by 



(J E (n,t)) =Tr[e iHt p e- im J E (n)] 



(4) 



Equation (Q| can be computed efficiently using the real- 
time 37 j finite-temperature j38| density matrix renormal- 
ization group [39|, [4(| (DMRG) algorithm introduced in 
PH . DMRG is essentially controlled by the so-called dis- 
carded weight e. We ensure that e is chosen small enough 
and that N is chosen large enough to obtain numerically- 
exact results in the thermodynamic limit. 

Non-equilibrium energy current — We start by study- 
ing a XXZ chain with two additional perturbations 
(dimerization A < 1 and a staggered field b > 0) which 
both render the system nonintegrable [1,[3(], 31 1. At time 
t = 0, two semi-infinite chains each being prepared in 
thermal equilibrium at temperatures Tj,^ are coupled 
by ho to an overall translationally-invariant chain. Ex- 
emplary results for {Je{ji, t)) are shown in Figure[T] The 
current at the interface n = saturates on a scale t ~ 1 
[note the definition of units via Equation ^)] irrespec- 
tive of the temperature difference Tl — Tr or the abso- 
lute values of Tl.r and regardless of the fact whether or 



3 




i i i i i i i i i i i i i i 

1 2 3 4 5 0.1 



Tr T r 

FIG. 3. (a) Temperature-dependence of the steady-state energy current (which becomes position-independent for t —¥ oo) 
of the integrable XXZ chain in the gapless (A = 0.5) and gapped (A = 2) regime. The right temperature Tr is varied at 
fixed Tl- (b) The curves at different Tl collapse if shifted vertically; thus, the T-dependence is of the simple functional form 
lmii^oo (Je(ti, t)) — /(Tl) — f (Tr). The linear thermal conductance ~ drf therefore completely determines the nonequilibrium 
current. A CFT approach [3^| predicts universal behavior f(T) = Cj^T 2 at low T which we verify (Insets) for the gapless XXZ 
chain (central charge c = 1) and the quantum Ising model (c = 1/2). The free fermion case A = can be solved analytically 
[2S| . |29| | at any T and provides a test for our DMRG numerics. 



not the system is gapped. The only exception is b > 
where ( Je(»i, t)) does not reach a finite steady-state value 
within the time scales accessible by our numerics [Figure 
QJc)] , again irrespective of the fact whether or not b opens 
a gap. We will now try to understand this in more detail. 

The time evolution of the local energy density h(n, t) = 
(h n {t)) of the XXZ chain (which for a homogeneous sys- 
tem might be a measure for an effective temperature) 
is shown in the Inset to Figure [lja). It does not reach 
a steady-state value but becomes increasingly smooth. 
This is not suprising since we are simulating a closed 
quantum system - but gives rise to the questions: (1) 
Why does the current saturate except for b > 0? (2) Are 
we actually modelling energy transport between reser- 
voirs whose temperature is fixed [27[? Both are reason- 
able if it does not matter over which length scale L the 
temperature difference Tl — Tr is applied; qualitatively, 
this should be the case if thermal transport properties of 
the chain are length- independent, i.e. if the thermal con- 
ductance G of a finite system does not decrease with its 
length L, or equivalently, if the conductivity a = GL of 
an infinite chain L — > 00 is infinite. More quantitatively, 
we conjecture an intimate relation between nonequilib- 
rium and linear response: The nonequilibrium energy 
current relaxes to a finite steady-state value if the linear- 
response thermal conductivity is infinite, i.e. if the Drudc 
weight D is nonzero. To support this we now proceed by 
calculating D. 

Linear response thermal Drude weight — The thermal 
Drudc weight is determined by the long-time behavior of 



the energy current correlation function 0, [i| 
Rc (J E (t)J E ) 



D = lim lim 

t— >oo N— >oc 



2NT 2 



Je 



J E {n) , (5) 



which can be readily computed using DMRG. Results 
are shown in Figure O For A = 1 and b = 0, Je is 
conserved, thus (Je^Je) = (Je(0)Je) [and the Drude 
weight can alternatively be obtained via Bethe ansatz 
see symbols in (a) and (b)]. The energy current 
correlation function of the dimerized chain seems to sat- 
urate at a finite value (or decays on a hidden large time 
scale). Our data thus indicates a nonzero Drude weight 
- which is interesting because the model is nonintegrable 
(previous results for D [1, 0, [l5| are controversial) . The 
deeper reason for this - the protection of the energy cur- 
rent by an unknown nonlocal conserved operator [3j - will 
be left as a subject for future work. Our data for b > 0, 
however, is consistent with D = [see Figure Hfc)] both 
in the gapless and gapped regime. Recalling that the 
nonequilibrium current relaxes to a nonzero steady-state 
value in all cases where 6 = 0, the observation of a finite 
(vanishing) Drude weight for A < 1 (b > 0) supports our 
above conjecture. 

Asymptotic current, homogeneous system — We now 
turn to study the temperature-dependence of the steady- 
state (position-independent) current. The result for the 
XXZ chain both in the gapless in gapped regime is illus- 
trated in Figure [3Ja) . The asymptotic current seem to 
be of a strikingly simple functional form: 



]im(J E (n,t)) = f(T L )-f(T R ) 

t— >OC 



(6) 



indicating a second intimate relation between nonequi- 
librium and linear response: The linear thermal conduc- 
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tance ~ drf(T) determines the steady-state noncqui- 
librium current at any Tl — Tr. Equation ([6]) can be 
established by varying Tr at fixed Tl; the correspond- 
ing curves collapse if shifted vertically [see Figure EJJb)] . 
The limiting behavior (both in the gapless and gapped 
regime) of f(T) is given by 



(7) 



Other details of / such as prcfactors or the crossover 
scale (which for A = and A 3> 1 is determined by 
the bandwidth or the size of the gap, respectively) in 
general depend on the model parameters. However, a 
recent conformal field theory approach [32I ] conjectures 
that the low-temperature behavior of a gapless system is 
universally given by f(T) = cf^T 2 with c being the CFT 
central charge, which follows intuitively from the version 
of the Stefan-Boltzmann law satisfied by a CFT [33| . We 
confirm this prediction for the XXZ chain (c = 1) as well 
as the quantum Ising model (c = 1/2); this is illustrated 
in the Insets to Figure GJfb). This is a highly nontrivial 
result because: (1) It is unclear why for a microscopic 
model whose equilibrium physics is governed by a certain 
low-energy field theory the very same field theory should 
describe the long-time behavior of the microscopic model 
in noncquilibrium (note that the behavior for < b < 
b c is not captured by the CFT!), and (2) Even linear- 
response transport properties (such as the Drude weight) 
are not determined by the low-energy theory alone but by 
a delicate interplay of conserved quantities protecting the 
current and dangerously irrelevant operators scattering 
the latter [H-i]. 

Our results for the dimerized chain at A = 0.8 are still 
consistent with Equation ([6]l . indicating that it might 
be a universal property of any translationally-invariant 
system. At smaller A and low T, we cannot reach time 
scales where oscillations of the current have died out com- 
pletely; we thus cannot exclude a different behavior of 
models which are strongly nonintegrable. 

The free fermion case A = can be solved exactly 
p8l [2^; Equation ([5]) reflects a noninteracting thermal 
Landaucr-Biittikcr formula. This analytic result can be 
used to test our DMRG numerics at any temperature [see 
the comparision in Figure [3^b) as well as in the Inset to 
Figure 0] . 

Inhomogeneous systems — We finally investigate sys- 
tems which are not translationally invariant. We start by 
studying the effects of different bulk parameters J n <o — 
J L , A„ <0 = A L , J„>o = J R , A„> = A R (and b„ = 0). 
If these parameters are chosen such that the renormal- 
ized Fermi velocities in the left and right halves coin- 
cide, backscattering due to the barrier (which is naturally 
present at the interface n = EH) can be tuned to zero 
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|4l| . and Luttinger liquid (LL) boundary physics 34 1 is 



absent. Results are displayed in FigureHJ indicating that 



FIG. 4. Steady state current for A = 1, b = but different 
bulk interactions Al = 0.3, A_r = 0.8 and different couplings 
Jl = 1, Jr = 0.808 to tune the backscattering at the interface 
to zero [ilj . The magnitude of the current is the same for 
T L = Ti,T R = T 2 and T L = T 2 ,T R = Ti. The curves at 
different T2 collapes, indicating that Equation ((6| might still 
hold. Inset: Linear thermal conductance of two identical XXZ 
chains connected by a strong barrier Jo = 0.002. The low-T 
behavior is consistent with r p 2 / K - 1 where K is the Luttinger 
liquid parameter [34|] . At A = we repoduce the exact result 
of [H for any T. 



Equation © might still hold. 

If two homogeneous XXZ chains J n =£o — 1,A„ = A 
arc connected through a barrier Jo, the linear thermal 
conductance G = Qtl {Je(^ — > 00)) \t l =t r =t is expected 
to feature a low-T power law T 2 l K ~ l with K being the 
LL parameter I34I For A = 0, Equation ((6]) still holds 
even if J < 1 [29|, and 2/K — 1 = 1 just reflects the 
asymptotic T 2 behavior of /. Our data for G in presence 
of interactions is consistent with y 2 /^- 1 . This is illus- 
trated in the Inset to Figure 2] where Jo is chosen small 
so that the scale on which LL boundary physics mani- 
fests becomes large (43|. Equation ([6]) still seems to hold 
above this scale. 

Summary - In this paper we provided evidence 
for an intimate connection between noncquilibrium and 
linear-response thermal transport properties of isolated 
infinite spin chains (or equivalently, interacting spinless 
fermions): (1) The energy current of a system which ini- 
tially features a temperature gradient Tl ^ Tr saturates 
to a finite value if the equilibrium thermal Drude weight 
D is finite, and (2) The value of the steady-state current 
at arbitrary Tl.r is completely determined by the lin- 
ear thermal conductance. Our data suggests that D > 
for a nonintegrable dimerized chain (or that the current 
correlation function decays on a hidden large time scale) . 
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